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Abstract 

The primary emphasis of this study has been to explain how modifying a cake 
recipe by changing either the dimensions of the cake or the amount of cake bat- 
ter alters the baking time. Restricting our consideration to the genoise, one of 
the basic cakes of classic French cuisine, we have obtained a semi-empirical for- 
mula for its baking time as a function of oven temperature, initial temperature 
of the cake batter, and dimensions of the unbaked cake. The formula, which is 
based on the Diffusion equation, has three adjustable parameters whose values 
are estimated from data obtained by baking genoises in cylindrical pans of var- 
ious diameters. The resulting formula for the baking time exhibits the scaling 
behavior typical of diffusion processes, i.e. the baking time is proportional to the 
(characteristic length scale) 2 of the cake. It also takes account of evaporation of 
moisture at the top surface of the cake, which appears to be a dominant factor 
affecting the baking time of a cake. In solving this problem we have obtained 
solutions of the Diffusion equation which are interpreted naturally and straightfor- 
wardly in the context of heat transfer; however, when interpreted in the context 
of the Schrodinger equation, they are somewhat peculiar. The solutions describe 
a system whose mass assumes different values in two different regions of space. 
Furthermore, the solutions exhibit characteristics similar to the evanescent modes 
associated with light waves propagating in a wave guide. When we consider the 
Schrodinger equation as a non-relativistic limit of the Klein-Gordon equation so 
that it includes a mass term, these are no longer solutions. 

1 Introduction 

It is an aphorism (at least to some physicists) that all natural phenomena can be under- 
stood and clarified in accord with general physical principles. The examples of such, 
that come immediately to mind are those related directly to the natural sciences. There 
are, however, numerous other examples, where the connection to physics is blatantly 
obvious and tacitly assumed, but is, nonetheless, usually ignored. One notable ex- 
ample is the culinary arts. In order to cook and do it well, in all likelihood, requires 
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little knowledge of the sciences in general and physics in particular. Indeed, most cook 
books are just that cook books: one merely follows the steps in a recipe, and one ob- 
tains the finished product. In general, no information is given detailing the physical 
principles on which a recipe is based, making it difficult, in many instances, to modify 
the recipe in non-trivial ways. 

Surprisingly, there appears to be a dearth of research related to the cake baking 
process. In fact, much of what appears in the literature is specifically related to bread 
baking. A notable exception is the work of Lostie and others 1 12|,| 1 1 1. In that study 
they have analyzed, in detail, the baking process of sponge batter during the first bak- 
ing period. 1 ' 2 As described in reference [ 12 1 the data used in their analysis have been 
obtained from highly controlled experiments in which the cake batter has been heated 
only from the top, with the sides and bottom of the pan being thermally insulated. 
Based on these experiments they have proposed a one dimensional model of the bak- 
ing process to predict the spatial dependence and time evolution of water content (water 
and vapor), temperature, gas phase pressure and porosity (the proportion of non-solid 
volume to total volume) of the cake batter. The proposed model has nine adjustable 
parameters whose values are determined by fitting the model to the data described in 
part I of their study. This study appears to be the most comprehensive in providing a 
complete understanding of the baking of sponge batter. 

The motivation for this article is to provide an answer, albeit a limited one, to the 
question of how the baking time of a cake varies when its recipe is modified, either by 
changing the dimensions of the cake or the amount of cake batter. In order to make 
the problem manageable we have restricted our analysis to one type of cake only, the 
genoise, one of the basic cakes of classic French cuisine. Although the technique for 
baking this cake is somewhat exacting, the ingredients should be readily available, 
facilitating attempts to reproduce our results. The theoretical basis of our solution is 
the Diffusion equation, and because of its intimate relation to Schrodinger equation, 
its solutions, in a certain sense, can be considered dual to solutions of the Schrodinger 
equation. We explore the dual interpretation of the solution to the cake baking problem 
and find that its interpretation from the perspective of the Schrodinger equation to be 
somewhat surprising and unsual. 

The outline of the article is as follows. First, we briefly discuss the quantitative 
aspects of cake baking. Next, we describe the experimental techniques and data used 
in the subsequent analysis. We, then, propose a simple model for the cake baking 
process and use it to obtain a semiempirical formula for estimating the baking time of 
the genoise. Finally, we interpret the model in the context of the Schrodinger equation 
and discuss the implications. 

1 The baking process involves several stages before the cake is completely baked. Their study only applies 
to the initial stage of the baking process. 

2 In reference [ 1 1 ] a comprehensive review of the existing research is given. 
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2 Cake Baking as a Diffusion Process 



2.1 Cake Baking from a Quantitative Point of View 

Before modeling the cake baking process, mathematically, we must first make more 
precise the imprecise measuring techniques that have, customarily, accompanied the 
instructions for baking a cake. A typical cake recipe lists the ingredients with amounts, 
a description of the technique used to prepare the cake batter, the baking time for 
the recipe, and usually some qualitative measure of determining whether the cake is 
baked. Until recently, amounts of ingredients have generally been given in units of 
liquid measure, like the cup or tablespoon. One can only surmise that the reason for 
this has been the scarcity of mass measuring scales in the home. The obvious reason 
for using dry measure (measuring by scale), rather than liquid measure, is to increase 
the likelihood of consistently producing a cake with the desired taste, texture, and 
appearance. Generally, the time necessary to bake a cake is accompanied by some 
subjective measure of determining whether the cake has baked to the proper degree of 
doneness. Such subjective measures include, when a toothpick thrust into the center of 
the cake comes out clean, when the cake shows a faint line of shrinkage from the sides 
of the pan, when the top springs back lightly, when you can smell the cake, etc. To a 
first approximation this criterion of doneness can be quantified by equating the degree 
of doneness to the temperature at the center of the cake, although one may imagine 
more complex criteria, which include other measurable quantities, like the temperature 
gradient. 

2.2 The Genoise 

The cake used in this study is the genoise, a basic cakes of classic French cuisine. 
The recipe 3 is adapted from that described in a book of Jacques Pepin 1 13 1. 4 As is 
typical, amounts of most ingredients are given in liquid measure. In Tabled we have 
converted these to dry measure. Qualitatively, we have found that the recipe prepared 
using traditional measuring techniques produces a cake consistent in texture, taste, and 
appearance with one prepared using dry measuring techniques. In fact, quantitatively, 
amounts measured and re-measured using traditional techniques differ, typically, less 
than two percent when compared to amounts measured by weight. According to the 
recipe the cake batter is placed in two cake pans, 8" in diameter by 1 \ deep, filling 
each pan | full. The cakes are then baked in an oven of temperature T a — 350° F for a 
duration of time between 22 and 25 minutes. In order to establish benchmarks for this 
study we have prepared the recipe, with exceptions as noted, filling two 8" cake pans to 
a depth of 1". Assessing the degree of doneness using traditional techniques we have 
observed the baking time to be, approximately, 17 minutes. The qualitative measure of 
doneness corresponds, quantitatively, to a temperature Tf (Tf = 203° F) at the center 

3 The only significant modification to the recipe is that the eggs have not been warmed over boiling water 
prior to mixing. 

4 Barham provides an excellent discussion about the physics and chemistry of baking sponge cakes, of 
which the genoise is an example 
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Traditional Measure Dry Measure (gm) 


eggs 


6 large 298 


sugar 


3/4 cup 176 


vanilla extract 


1/2 tsp 2 


all purpose flour 


1 cup 144 


butter 


3/4 stick 114 



Table 1 : The proportion of ingredients for the genoise is given in traditional measure 
as well as dry measure. 



of the cake. The recipe has been prepared several times with baking times varying by 
no more than two minutes. 

2.3 Theory from a Naive Perspective 

Initially, we were believing that we could model the cake baking process, naively, as a 
simple diffusion process. 5 Although the model is grossly inadequate, it has provided 
the basis of the final model and is, thus, presented. This model derives from a solution 
of the Diffusion equation |9 1, 



at 

where T = T(t, r, 8, z) is the temperature of the cake at time, t, at location, (r, 8, z) 
(expressed in cylindrical coordinates), within the cake, and T> is the diffusion coef- 
ficient (assumed constant) for the cake batter. We assume the cake batter to be in a 
cylindrical pan of radius, R, and thickness, Z, at initial temperature T c , and baked in 
an oven of constant temperature, T . Assuming these initial conditions we solve the 
Diffusion equation for the temperature, T, which because of azimuthal symmetry is 
independent of 8, 

T(t,r,z)=T + (T c -T )^l V [1 " M) "U m » , (2.2) 

7T f— ' nx m fi 
n— l,m— 1 

where x m ,o is the m-th root of the zeroth order Bessel function, and the functions $ mn 
are the normalized eigenfunctions of the Diffusion Equation, 

'> nnz^ Jo( x "j. or ) 



dT 



$m,n = — — —j= Sin(— -) 



RJ\[x m fi)\fZ Z x mt oJi(x m fi) (2 3) 

xex P [-{(^) 2 + (^) 2 Pi 

In Equation 12.21 the first term to the right of the equality is the steady state solution, 
i.e. a solution of the Laplace equation. The second term is the solution particular to 
the boundary conditions, i.e. vanishing at the surface bounding the cylinder. Only the 

5 Klamkin has considered such a model for estimating the scaling behavior in the cooking time of a 
roast H51 . 
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Bessel functions of zeroth order appear because of azimuthal symmetry. The other 
factors are required to reproduce the initial temperature on the interior of the cake at 
t = 0. 

Typically, the initial temperature of the cake batter T c w 80° F, corresponding 
to the temperature at which the cake batter is prepared. Thus, the initial and final 
temperatures of the cake and temperature of the oven satisfy the relationship 

T/ ~ Tc < 1 . (2.4) 
T -T c ~ 

For temperatures, Tf which satisfy Eg. 12.41 the infinite series, Eg 12.21 evaluated at the 
center of the cake, can be approximated by a single term 

r(t/,0, Z/2) =T f = T + (T c - T )F , (2.5) 

where 



F 



^ ^-{({I'+mpt/i 

if Ci < 1 and C 2 < 1 

|ex P H(§) 2 }^/] 
if Ci > 1 



-^eM-{(( X -je) 2 )}Vt f ] 
if C 2 > 1 



Here 



C x = 



^-eM-{C-^) 2 )}Vt f ] 



C 2 = lexph{(|) 2 }^ / ], 

7T Z 

and tf is the amount of time for the temperature of the center of the cake to increase 
to temperature Tf. This simplification to Eg. 12.21 is most readily verified by direct 
calculation. 

We can solve Eguation l2.5l explicitlv for the baking time, tf, 



(2.6) 



(2.7) 

ED 



x l,0 I 7T 



— oc < 



™ i o 



if C x < 1 and C 2 < 1 



if Ci > 1 



if C 2 > 1 



(2.8) 
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where 

x\ n = the first root of the zeroth order Bessel Function 

(2.9) 

w 2.40 



The constant of proportionality depends only on the unknown diffusion constant, T>, 
which can, in principle, be estimated from the baking time of a given recipe. In any 
case, if the baking time is known for a cake of specified dimensions, then Eg. 12. 81 c an be 
used to estimate the baking time of a cake of arbitrary dimensions. For practical appli- 
cations it is useful to express Eq. 12. 8 1 with its numerical constants appearing explicitly, 
in which case we have 

2.34 , J_ 
D 2 ' Z 2 

if d < 1 and C 2 < 1 

l 

ifCi>l , (2.10) 

2.34 

if C 2 > 1 



where D = the diameter of the cake , 6 ' 7 

We use Ea. l2.10l to estimate the baking time of a cake which is 2.1" radius and 4.0" 
depth. The calculation assumes the same initial conditions as apply in the preparation 
of the basic recipe described previously. 8 The predicted baking time is approximately 
85 minutes. The actual baking time, based on an internal temperature of 203° F, is 26 
minutes. 



oc < 



3 Experiment 
3.1 Procedure 

In order to understand why the previously described model of the cake baking process 
has failed, we have collected a variety of data for different amounts of cake batter baked 
in cylindrical cake pans of various dimensions. Seven cakes have been baked whose 
dimensions range from 2.1" to 6.5" radius and from 1.0" to 4.0" depth and whose 

6 It is interesting to note that if the cake is rectangular, the term 2.34/D 2 is replaced by 1/X' 2 + 1/Y 2 , 
where X and Y are the length and width of the cake pan. 

7 The reader should be aware that some formulae are expressed in terms of diameter rather than radius. 
The reason is to make transparent the relationship between these formula those applicable to rectangular 
cake pans. 

8 It is worth mentioning that standard cake pans are not readily available with these dimensions. This has 
necessitated substituting a coffee can, filled with batter to a depth of 4.0", for a cake pan. In fact, genoises 
are commonly not baked greater than 2" in depth. 
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baking times range from 17 to 42 minutes. The data consist of the mass of each cake 
before and after baking, the diameter and depth of the cake, the depth of the cake after 
baking, and the temperature at the center of the cake recorded at one minute intervals 
until the cake is baked. The experimental setup for recording the internal temperature 
of the cake is shown in Figure^ In most ovens the temperature of the oven is measured 
by a temperature probe which is close to the oven walls, and consequently, relatively 
distant from whatever is being baked in the oven. The temperature at the oven walls is 
typically higher than the temperature at the surface of what is being baked. In this study 
temperatures have been monitored next to the cake pans, with the heating element of the 
oven regulated manually to maintain the desired oven temperature. This has been done 
to minimize convective effects so that transfer of heat to the cake can be approximated 
as a process involving only conduction. 

The experimental procedure involves the following: 

1. measure the mass and diameter of the cake pan; 

2. prepare the cake batter; 

3. fill the cake pan to the desired depth; 

4. measure the mass of the cake batter; 

5. place the cake in the preheated oven and record its internal temperature at one 
minute intervals until the cake is baked, i.e. until the cake's internal temperature 
reaches 203°. 9 

6. after the cake is baked, re-measure its depth and mass. 

It is in order to make some comments about the experiment. The baking time of 
individual cakes is extremely sensitive to the depth of the cake batter. The reason is 
that the batter is relatively shallow, as is typical when baking cakes of this type. For 
example, if the cake batter is approximately 1 ~ deep, an error in measurement of | 
is a 12% error. Even when using extreme care to make depth measurements at the cen- 
ter of the cake it is difficult to insure that the depth remains constant across the surface 
of the cake. This problem is exacerbated in the cakes of larger diameter used in this 
study. The fact that cakes rise while baking, in itself, is not a problem, but identical 
cake recipes baked on different occasions do not rise precisely the same way. These un- 
certainties help to clarify why for many cake recipes a range of baking times are given 
whose extremes differ by 15% or more. It seems reasonable to assume that this type 
of random behavior is more exaggerated in cakes of larger radius resulting in greater 
uncertainties in their baking times. There is another practical consideration which is a 
potential source of error in the experiment. Because of the shallow depth of the batter 
relative to the size of the temperature probe it is difficult to insure that the probe is in- 
serted accurately into the center of the cake. On the other hand, if extremely controlled 
experiments are needed for making accurate theoretical predictions, then from a prac- 
tical perspective the results would be of little use. One final comment, as it has been 

'Temperature measuring instruments are located outside of the oven, so that temperature data can be 
collected without opening the oven door. 
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impractical to automate the data collection process, specifically the temperature/time 
measurements, the actual data collection has been tedious. 




Figure 1: Shown are the oven, the baking pan, and two thermocouple sensors and 
temperature measuring instruments. One sensor monitors the internal temperature of 
the cake, and the other the temperature of the oven. 

In Figure I2 we show a plot of temperature as a function of time at the center of a 
cake 4.0" radius and 1.0" depth. In addition, we also show the theoretical fit to the data 
based on Eg. 12.21 The summations in Ea. l2.2l have been truncated at n, to = 17. 10 It 
appears that the model based on the Diffusion Equation is inadequate in accounting for 
the features of the data. In Figure[3]we show a plot of temperature as a function of time 
of a cake 2.1" radius and 4.0" depth. In this case, also, the theoretical curve appears to 
fit the data inadequately, suggesting that a drastically different theoretical model from 
that proposed, initially, may be required to estimate baking times. 

3.2 Revised Model 

After analyzing the data both qualitatively and quantitatively, we have made several 
observations and drawn a number of conclusions. We have observed that for a baked 
cake, not only does the volume of the cake increase, but also the mass of the cake de- 
creases. 11 It also seems apparent that the model based on Eg. 12. 21 exclusively, provides 
a poor fit to the data of temperature verses time for the cakes of various dimensions. 
With hindsight it has been naive to believe that a model based, only, on Eq. 12.21 can 
account for the features of the data. The cake batter is a viscous liquid at the start of 
the baking process and a solid at the end. Thus, not only conduction, but also convec- 
tion, and even radiation contribute to the baking process, each to varying degrees. In 

10 It has been necessary to truncate the sum at such large values of m and n to obtain reasonable estimates 
of temperatures close to the initial temperature of the cake. 

"We have observed that typically the mass of a baked cake is less by approximately 10%. 
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Cake Dimensions - 4.0 inch radius, 1 .0 inch depth 
220 | 1 1 1 1 1 1 1 




80 I ^ 1 1 1 1 1 1 1 1 

2 4 6 8 10 12 14 16 18 

Time (min) 

Figure 2: The temperature of a genoise 4.0" radius and 1.0" depth recorded at one 
minute intervals is shown along with the theoretical fit based on Eg. 12. 21 

addition, the parameters associated with these methods of heat transfer, thermal con- 
ductivity, density, specific heat, the diffusion coefficient, and convection coefficient, to 
some extent, depend on temperature. In any case, assuming that the diffusion coeffi- 
cient is constant, as in Eg 12.21 is a simplistic assumption. To improve the model we 
make the following assumptions and revise the model accordingly: 

1 . conduction is the primary mechanism of heat transfer; 

2. the thermal conductivity of the cake batter varies during the baking process, pri- 
marily, because of moisture evaporation at the top, cake surface. 

In the revised model moisture evaporates at the top surface of the cake causing the 
cake batter to dry out so that, consequently, the diffusion coefficient decreases in value. 
The cake batter consists of two homogeneous mixtures, each with its own diffusion co- 
efficient. In Figure|4]we diagrammatically represent the process. As baking occurs the 
interface between the two homogeneous mixtures of cake batter (dry on top and moist 
on the bottom) moves from the top to the bottom of the pan. The quantity a denotes the 
interface between the two mixtures at some arbitrary time. Here, T> and V are the dif- 
fusion constants of the dry and moist mixtures, respectively. In order to reduce the time 
for numerical computations we have implemented the following simplification in ap- 
plying the model in practice. We have assumed that the cake batter remains uniformly 
moist during the time interval < t < t\ (Interval I). Then, during the time interval 
ii < t < ti (Interval II) the cake comprises two homogenous mixtures, dry batter in 
the top half and moist batter in the bottom half of the pan, i.e. the interface remains 
fixed at a = Z/2. At time t-i and until it completes baking at time tf, i.e. t2 < t < tf 
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Figure 3: The temperature of a genoise 2.1" radius and 4.0" depth recorded at one 
minute intervals along with the theoretical fit based on Eg. 12. 21 

(Interval III), the cake is uniformly dry. 12 

For time t < t\ and ti < t < if the time evolution of the temperature is calculated 
using the eigenfunctions of the Diffusion equation given in Eq. 12.21 substituting the 
diffusion coefficient V or V , as appropriate. For time t,t\ < t < i2, however, we re- 
quire eigenfunctions obtained from solving the Diffusion equation for two contiguous, 
homogenous media whose interface is located at z = a. 

The unnormalized eigenfunctions, ^ m j, of eigenvalue, E m j, are given by 



= exp(-E mj t) x J ( ^I^L j Z mj (z) . (3.1) 



Here 

Z_,fVl = J 

B mj sin(A^ [Z - z)) a < z < Z 



-7 I \ - J Am J sin ( A »u z )( z ) < z < a 

^rnj\ Z ) — \ n try \\ ^ r, 1 



The coefficients A m j and £? m:) are determined by the normalization condition and the 
boundary conditions at z = a. The quantities X m j and \' m j are related by 



(^) + • (3-3) 



and satisfy the eigenvalue equation 

tan(A mj a) _ ta,n(X' mj (Z - a)) 



(3.4) 



12 Any time dependence of the diffusion coefficient is assumed to result from that of the thermal conduc- 
tivity, only. Thus, we neglect possible time dependence in the mass density or specific heat. This simplifies 
the analysis and avoids physically unrealistic, temporal discontinuities in the temperature which result from 
discontinuous changes in the diffusion coefficient. 
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Moisture Flow 



Figure 4: As heat flows into the cake, moisture evaporates through the top surface 
of the cake resulting in a two component system, one with a diffusion coefficient V, 
corresponding to the moist batter in the bottom portion of the pan and another with 
diffusion coefficient V, corresponding to the dry batter in the top portion of the pan. 



The eigenvalue equation is obtained from the boundary conditions at the interface 
separating the two homogeneous mixtures of cake batter: 



Z(z)\ z=a - =Z'{z)\ t=a+ , 



(3.5) 



and 



V 



dZ(z) 



dz 



= V 



,dZ'{z) 



dz 



(3.6) 



The latter equation is the result of applying the continuity equation at the interface. 
Here T> > V because the moist cake batter thermally conducts better than the dry cake 
batter. Expressed as a function of A' ■, 



\2 

mj 



V 

— A' 2 

V A rnj 



V-V 
V 



*£m,0 



(3.7) 



The solutions of the eigenvalue equation exhibit the following properties: 



1 A 2 < A' 2 • 

" raj — mj ' 



2. there are no solutions for A^ < ; 

3. the derivative of an eigenfunction is discontinuous at z = a, 

4. Some eigenfunctions are oscillatory in the region < z < a and decaying in the 
region a < z < Z. 



Some of these properties have interesting consequences, especially when interpreted 
in the context of the Schrodinger equation. We defer further discussion of these until 
Section|5] 

We now prescribe the recipe for implementing the revised model in practice. 
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1. Specify values for the adjustable parameters T>, V, t\, and t 2 . 

2. Substitute the temperature of the oven, T a , the initial temperature, T c of the cake, 
the radius, R, of the cake pan, and the depth, Z, of the cake as well as the value 
of V into Eg. 12. 21 This gives the internal temperature of the cake, T] (t) for time 
t, < t < t x . 

3. Tj(ii) becomes the initial temperature of the temperature distribution for times 
t, ti < t < t 2 . 

4. Expand Tj(ti) in terms of the eigenfunctions, f mj -, Eg. 13. II to obtain the tem- 
perature distribution Tjj(t) during the time interval t\ < t < t 2 . 

5. Tufa) becomes the initial temperature of the temperature distribution for times 
t, h < t < tf. 

6. Expand Tjj(t 2 ) in terms of the eigenfunctions, $ mn , Eg. 12.31 to obtain the tem- 
perature distribution Tjjj(t) during the time interval t 2 < t < tf. The diffusion 
coefficient V is substituted for T> in the $ mrl . 

In practice, the evaluation of the infinite series for the temperature distributions in In- 
tervals I,H, and III requires that each series be truncated at a finite number of terms, 
depending on the degree of accuracy desired. For Interval I, assume that Eq. 12.21 is 
terminated for integers M, (m < M) and N, (n < N). In order to achieve approx- 
imately the same degree of accuracy in Interval III, we apply the same values to the 
corresponding infinite series. For Interval II, the temperature distribution is expanded 
using all & m j, also, for values of m < M, but selection of the truncation for j is more 
complicated. All values of j such that 

Am , < mirni (3 . 8) 

Z — a 

are used. We note that solving for all eigenvalues A mj is nontrivial because there is 
not a one to one correspondence between those in Interval II and those in Intervals I 
and III. 

We have used these expansions to model the temperature at the center of the cake 
for four representative, sets of data. We summarize the results in Tableland Figures|5] 
through[8] 13 The data have been selected to span the range of cake dimensions. As is 
apparent in Figures |5] through [8] it is possible to obtain a reasonable fit to each set of 
data by varying the four adjustable parameters in the model, T>, V , t\,t 2 . The primary 
purpose of constructing these models has been to provide incite into the baking process 
so that quantitative estimates can be obtained for the baking time of cakes of various 
dimensions. These results are consistent with the qualitative aspects of the revised 
model. It appears, in accordance with Table [2] that cakes of larger depth, although 
taking longer to bake, have an effective diffusion coefficient which is greater than cakes 

13 Computations have been performed using Maple software running under Windows XP on a one GHz, 
dual processor, Pentium III computer. The optimization of adjustable parameters in the model has been 
achieved by trial and error, as there has been no simple or straightforward way to interface the model calcu- 
lations to optimization software. 
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Data 


Adjustable Parameters 


Accuracy 


Radius (in) 


Depth (in) 


tf (min) 


P(in^/min) 


X>'(in^/min) 


ti (min) 


t 2 (min) 


M N 


4.0 


1.0 


17 


.029 


.0035 


2 


2 


10 10 


2.1 


4.0 


26 


.012 


.047 


10 


oo 


10 10 


2.1 


2.0 


20 


.025 


.010 


1 


15 


10 10 


5.0 


1.8 


41 


.035 


.005 


2 


2 


10 10 


6.5 


1.6 


35 












5.0 


1.0 


20 












2.1 


1.8 


19 













Table 2: Shown is a summary of the data used in this study. A representative sample 
of these data has been selected for purposes of modeling. The values M and N have 
been selected to insure that the theoretical curves are relatively smooth. The adjustable 
parameters have been determined by trial and error. 



of smaller depth. On the other hand, cakes of smaller radius seem to have an effective 
diffusion coefficient which is greater than cakes of larger radius. 

In cakes of less depth moisture can evaporate from the top surface of the cake 
relatively quickly causing the cake to dry out. This causes a decrease in the diffusion 
coefficient because of lack of moisture content. For an amount of baking time, t, we 
hypothesize that the moisture content within the cake scales, approximately, as s tz , 

st* = • (3.9) 

This type of behavior is typical in random walk processes since diffusion times scale 
as some ( characteristic length) 2 . 

Moisture, also, diffuses from the sides of the cake, inward. In situations when 
the radius of the cake is relatively less than its depth, moisture in the interior of the 
cake may remain constant, or possibly increase, even though moisture is continuing 
to evaporate from the top cake surface. For an amount of baking time, t, we further 
hypothesize that the moisture content scales, approximately, as s rz , 

< 31 °> 

We acknowledge that these hypotheses regarding the scaling behavior of the mois- 
ture content are highly speculative, and that the degree to which they can be justified 
rests on how successfully they can be applied to predicting the baking time of cakes of 
various dimensions. 



4 Estimating the Baking Time of a Cake 

We now utilize the revised model of Section l3~2l to obtain an explicit formula for the 
baking time of a cake. 14 In the revised model we have assumed that as a cake bakes 

14 In this Section the diameter, rather than the radius, is used in all formulae. 
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Cake Dimensions - 4.0 inch radius, 1 .0 inch depth 
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Figure 5: Shown is a theoretical fit to temperature as a function of time at the center of 
a cake, 4.0" radius and 1.0" depth. The relevant parameters are given in Table|2] 

its diffusion coefficient changes, primarily, because of evaporation of moisture at the 
top surface of the cake; however, for a given cake we can calculate, using Eq. 12.51 
an effective diffusion coefficient, 2? e ff, which is constant during the baking process. 
Furthermore, in accord with the revised model the functional form of 2? e ff should be 
consistent with Eas. 13.91 and 13.101 If the functional form is known, it provides a re- 
lationship between the baking time and physical dimensions of the cake, making it 
possible to solve for the baking time as a function of the cake dimensions. 

For each of the seven data entries in Table |2j we estimate an effective diffusion 
coefficient using Eq. 12.51 We, then, fit to these data an expression for 2? e ff of the 
following form, 



The adjustable parameters are estimated to be ao = 2.57 x 10~ 2 , a\ = 2.07 x 10~ 2 , 
and et2 = 19.2. The functional form of Ea. l4.1l has been chosen primarly because it 
can be interpreted in a relatively straightforward and physically intuitive manner. We 
have tried other formulae with more physical appeal, i.e. different dependencies on 
the baking time and cake dimensions; however, when inverted to estimate baking times 
they yield multiple roots for a cake of specified dimensions. Generally, these additional 
roots have no obvious physical interpretation. We have also varied the functional form 
of Eq. 14. II among those which are not plagued with multiple roots for the baking time 
and have found the calculated baking times not to be drastically different from those 
obtained using Ea.O 

The effective diffusion coefficient depends explicitly on the baking time, the di- 
ameter of the pan, and depth of the cake batter. In Figures |9] and ^0] we show the 



f eff = o-a ~~ a i ex P — 0-2 — 
L W 




(4.1) 
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Cake Dimensions - 2.1 inch radius, 4.0 inch depth 
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5 10 15 20 25 30 
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Figure 6: Shown is a theoretical fit to temperature as a function of time at the center of 
a cake, 2.1" radius and 4.0" depth. The relevant parameters are given in Table|2] 



Cake Dimensions - 2.1 inch radius, 2.0 inch depth 
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Figure 7: Shown is a theoretical fit to temperature as a function of time at the center of 
a cake, 2.1" inch radius and 2.0" depth. The relevant parameters are given in TableEl 
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Figure 8: Shown is a theoretical fit to temperature as a function of time at the center 
of a cake, 5.0" radius and 1.8" depth. The relevant parameters are given in Table|2] 



experimental values of 2? e ff and the theoretical fit. Rather than plotting them as a func- 
tion of the baking time, diameter of the pan, and depth of the batter, we have plotted 
them as a function tjjZ 2 and D 2 jZ 2 . The adjustable parameters in Eg. 14. II have a 
simple, physical interpretation. If values of tf, D, or Z are such that the exponential 
becomes vanishingly small, the cake batter should be maximally moist and 2? e ff is ao. 
On the other hand, if their values are such that the exponential approaches one, the 
cake should be totally dry with the effective diffusion coefficient being ao — a\. The 
parameter a-i relates to how readily moisture evaporates from the cake. These three 
parameters should also depend on the initial moisture content of the cake batter. 

Combining Eqs. l2.5l and l4.ll we obtain an implicit formula for the baking time as 
a function of depth, diameter, oven temperature, initial temperature of the cake batter, 
and final temperature of the cake: 



' 4 <0 , n 2 



K = -V eB t f < 



if Ci < 1 and C 2 < 1 
if Ci > 1 

ifC 2 > 1 



(4.2) 
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Figure 9: Shown are the experimental values of the effective diffusion coefficient and 
the theoretical fit to these values. The independent variable has been selected based on 
the discussion of Section lT2l 
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Figure 10: Shown are the experimental values of the effective diffusion coefficient and 
the theoretical fit to these values. The independent variable has been selected based on 
the discussion of Section l3T2l 
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Baking Time as a Function of Depth and Diameter 




Figure 11: The theoretical and measured baking times of cakes of various depths and 
diameters are shown. 



where K is obtained from 



The conditions C\ and C2 are defined by Eqs. l2.7l and l2.7l . 

We use Eas. 14.21 and !4.3l to estimate baking times of cakes of various dimensions 
used in this study and compare them to the actual baking times. These results are 
presented in Figure[^ an d Tabled It should be noted that Eqs. l4.2l and l4.3l exhibit the 
scaling behavior that is typical of diffusion processes, i.e. baking times are proportional 
to (characteristic length) 2 . 

Philosophically speaking, a semi-empirical model should accomplish more than 
regurgitate what is already known: 

1. it should have predictive ability; 

2. it should also be valid in that it is consistent with established physical principles; 

3. it should not have an inordinately large number of freely adjustable parameters. 



7T Zl.u JUgl^Oj 

8 

if Ci < 1 and C 2 < 1 
if d > 1 



ci,o Ji(zi,o) 
2 

if C 2 > 1 



(4.3) 
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if (min) 


ttheor (min) 


flinreg (min) 


Diameter (in) 


Depth (in) 


17 


16 


21 


8.0 


1.0 


26 


26 


29 


4.1 


4.0 


20 


20 


18 


4.1 


2.0 


41 


45 


30 


9.9 


1.8 


35 


41 


37 


13.0 


1.6 


20 


16 


25 


9.9 


1.0 


19 


19 


17 


4.1 


1.8 



Table 3: Baking times and those predicted by the theoretical model are presented for 
the diameters and the depths of cakes used in this study. For comparison the baking 
time predicted by linear regression is also presented. 

Consider the four hypothetical cakes whose dimensions and predicted baking times are 
given in Table|4]For various reasons these cakes are impractical to bake. Nonetheless, 
they do provide an arena for qualitatively exploring the theoretical basis of Eqs. 14.21 
and !4.3l Even to someone with little cooking experience it is not surprising that cake 
number two requires 17 minutes to bake, given that a cake 8" diameter and 1" depth 
requires 17 minutes baking time (See Table[3}. It may be surprising that cake number 
three requires so little time to bake, given how long it takes cake number two to bake. 
The reason is that the moisture content of cake number three remains high because 
very little evaporation takes place during the cooking process. As a result, the diffu- 
sion coefficient remains relatively large so that the cake bakes quickly. Cake number 
four requires nearly 24 hours to bake. It is interesting to note that the baking time of 
cake number four reflects the scaling property of diffusion processes, i.e. its physical 
dimensions are thirty times those of cake number one and is therefore predicted to 
require 900 times as long to bake. 

The semi-empirical model for the baking time has three adjustable parameters. 
From a practical perspective one may ask if performing a linear regression, which also 
has three adjustable parameters, would produce a model with as good or possibly better 
predictive ability than the model presented. The result of a linearly regressing baking 
time on diameter and depth is 

t f = -3.00 + 2.34D + 5.74Z. (4.4) 

In Table |3 we present the baking times predicted by Eq. 14.41 for comparison to those 
predicted according to theory. 

Although, strictly speaking, our analysis has been applied to the genoise we can 
apply the results, qualitatively, to baking in general. The difference between cake bat- 
ter and pastry dough is not so much a matter of different ingredients, but rather the 
proportions of the same ingredients. Both typically are composed of flour, fats, and in- 
gredients which contain water, such as eggs, milk, juices, etc. Typically, pastry dough 
contains a very small amount of liquid in proportion to other ingredients. This explains 
why it takes such a long time, relatively speaking , to pre-bake a pastry shell (15 to 30 
minutes, for example), even though, pastry shells are generally, not much more than | 
to i in depth. Cookie dough which lies between pastry dough and cake batter, with 
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Cake 


ttheor (mill) 


Diameter (in) 


Depth (in) 


1 


1.56 


1.00 


1.00 


2 


17.0 


30.0 


1.00 


3 


1.82 


1.00 


30.0 


4 


1.40 x 10 3 


30.0 


30.0 



Table 4: The predicted baking times of four hypothetical cakes are presented. 



respect to moisture content, also requires a relatively longer time to bake than cake 
batter. Cookies require typically between ten and fifteen minutes to bake, even though 
they are usually not more than ^ thick. 

5 An Irrelevant but Intriguing Digression 

The similarity between the Diffusion equation and the Schrodinger equation is strik- 
ing. We replace t in the Diffusion equation by it /h to obtain the Schrodinger equation. 
The duality between solutions of the Diffusion equation and Schrodinger equation is 
typically not insightful since the corresponding solutions of the Schrodinger equation 
are those in the absence of a potential. Consider, however, the eigenfunctions, Eg. 13. II 
They are dual to the eigenfunctions of a particle confined to a cylindrical box. What is 
atypical is the mass of the particle assumes different values in the regions < z < a 
and a < z < Z. The mass, /i, is related to the diffusion coefficient, T>, /i = ft 2 / (22?). 
It is common for situations to exist where the diffusion coefficient varies with position; 
however, it is less common for the mass of a particle to vary with position. Nonetheless, 
there do exist situations where the mass of a particle, or quantum mechanical system, 
may be different in different regions of space. Quantum systems with spatially depen- 
dent effective mass have been used, extensively, in modeling the physical properties 
of various microstructures and semiconductor interfaces in condensed matter. 1 1] Such 
specialized models have been applied to studying the electronic properties of semi- 
conductors |4|, quantum wells and quantum dots [8|,[14|, 3 He clusters |3|, quantum 
clusters |5|, and graded alloys and heterostructures. |6|,[7|,| 15|. 

If one is guided by the seeming duality between the Diffusion equation and Schrodinger 
equation, then the Schrodinger equation is 

This insures that the Schrodinger equation conserves probability so that the continuity 
equation remains true, 

dP 

V.J = --, (5.2) 
where the current density, J, is given by 

J = ^[V>*VV>-(VV>*)V>] , (5.3) 
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and the probability density, P, is given by 



P = ip*ip . (5.4) 

If fi is placed to the left or right of the derivatives, except for specific combinations 
of such, there will be sources/sinks of probability, as can be shown by direct calcu- 
lation. A consequence which results from the unusual placement of /i in Eq. 15.11 is 
that when the mass, /i, discontinuously changes values between two regions, the spa- 
tial derivative of the wavefunction is discontinuous at the interface between the two 
regions. This effect is apparent in Eq. 13.51 whenever T> ^ V . This effect is demon- 
strated in Figures E]^] and^J which are examples of the z -dependent part of the 
eigenfunctions of the Diffusion equation given in Eq. 13.21 

Another apparent effect which is evident from a detailed analysis of Ea. l3.7l is that 
for sufficiently large values of m and sufficiently small values of j there are solutions 
with y? m j < 0. The implication is that some of the eigenfunctions are hybrid in that 
for < z < a they decay and for a < z < Z they oscillate. The corresponding 
eigenfunctions of the Schrodinger equation also exhibit this same behavior, decaying 
in the region of smaller mass and oscillating in the region of larger mass. This effect 
is shown in FiguresE]^] and 1 141 In Figure^Jthe eigenfunction, whose value of j 
is the largest of three examples, is oscillatory in both regions. This effect can also be 
shown to persist under more general boundary conditions than those imposed by the 
cake pan. 15 From the quantum mechanical perspective the values of m are associated 
with momentum tangent to the interface separating the regions of different values of 
mass. Thus, if the momentum tangent to the interface is sufficiently large, the wave 
function damps out in the region of less mass. This is reminiscent of the evanescent 
modes associated when light propagates along a wave guide. 

If one considers the Schrodinger equation as the non-relativstic limit of the Klein- 
Gordon equation, then one should also include the term mc 2 on the left side of Ea. lBTl 
It can, then, be shown that the hybrid eigenfunctions are no longer solutions. 

6 Conclusions 

The primary emphasis of this study has been to explain how modifying a cake recipe 
by changing either the dimensions of the cake or the amount of cake batter alters the 
baking time. Our analysis has been restricted to one particular cake, the genoise. Based 
on data collected and ensuing analysis we have concluded that conduction is the pri- 
mary mechanism of heat transfer and that the Diffusion equation provides a theoretical 
framework adequate for describing the baking process; however, the diffusion coeffi- 
cient does not remain constant during the baking process. From detailed analysis of 
the data we have concluded that the diffusion coefficient changes during the baking 
process, principally, because of evaporation of moisture at the top surface of the cake. 
Consequently, we have approximated the cake baking process as one in which the dif- 
fusion coefficient, effectively, assumes a constant value during the baking process. Its 

I5 The effect seems only to be a consequence of the fact that the mass assumes different values in two 
contiguous spatial regions. 
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The Normalized Eigenfunction Z%^\(z) 
1 1 — — T 




0.2 0.4 0.6 0.8 1 



z 

Figure 12: The eigenfunction Ze.i exhibits both decaying and oscillatory behavior in 
regions < z < 1/2 and 1/2 < z < 1, respectively. In addition, ^ is discontinu- 
ous at z = a, where a = 1/2. 

value depends on the diameter, depth, and baking time of the cake, all of which affect 
the moisture content of the cake. We, then, have proposed a semi-empirical formula 
for the effective diffusion coefficient consistent with the physical principles underlying 
the model. The formula has three adjustable parameters whose values have been esti- 
mated from the data. We have inverted the formula to obtain the baking time of a cake 
as a function of diameter and depth. The resulting formula exhibits scaling behavior 
typical of diffusion processes: the baking time scales as the (the characteristic size ) 2 
of the cake. Although we have not solved the problem of how long is required to bake 
a cake, in general, we have offered a qualitative explanation of those factors which 
most importantly determine baking times. In addition to the dimensions of the cake 
we have proposed that the moisture content of a cake is the dominant factor affecting 
its baking time. This proposition explains why pastry dough and cookie dough require 
such a relatively long time to bake in comparison to cake batter, after accounting for 
differences in their physical dimensions. 

As part of the analysis and modeling process we have required solutions of the Dif- 
fusion equation for two cylindrical media of the same diameter, separated by a com- 
mon circular boundary. The diffusion coefficient assumes a different constant value 
in each medium. These solutions are interpreted naturally and straightforwardly in 
the context of the Diffusion equation; however, when interpreted in the context of the 
Schrodinger equation, they are somewhat peculiar. The solutions describe a system 
whose mass assumes different values in two different regions of space. Furthermore, 
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The Normalized Eigenfunction Zq^{z) 
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Figure 13: The eigenfunction, exhibits both decaying and oscillatory behavior in 
regions < z < 1/2 and 1/2 < z < 1, respectively. In addition, d e z 2 is discontinu- 
ous atz = a, where a = 1/2. 

the solutions exhibit characteristics similar to the evanescent modes associated with 
light waves propagating in a wave guide; however, when we consider the Schrodinger 
equation as a non-relativistic limit of the Klein-Gordon equation so that it includes a 
mass term, these are no longer solutions. 
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The Normalized Eigenfunction Z 6 ^(z) 
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Figure 14: The eigenfunction, Ze.3 exhibits only oscillatory behavior in the region 
< z < 1. In addition, d 6 z 3 is discontinuous at z = a, where a = 1/2. 
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